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ABSTRACT 


A  measurement  technique  has  been  developed  which  makes  possible  the 
determination  of  the  acoustic  transmission  loss  for  various  paths 
between  a  source  and  a  receiver  at  relatively  short  ranges  in  the 
ocean.  The  method  is  based  upon  a  technique'  called  Time  Delay 
Spectrometry  (TDS).  The  theory  for  utilizing  TDS  to  perform 
measurements  on  an  underwater  test  range  with  a  homogeneous  water 
column  is  presented.  The  experimental  parameters  concerning  test 
gecmetry  together  with  criteria  for  the  source  signal  are  discussed  and 
the  limitations  are  described.  A  computer  program  written  in  HP  BASIC 
is  the  tool  in  a  multipath  measurement  system  where  emalog  sound 
signals  recorded  on  tape  are  processed  by  am  HP  3561A  Dynamic  Signal 
Analyzer  under  program  control. 
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I.  INTROPaCTION 

A.  BACKGROUND 

The  Naval  Undersea  Warfare  Engineering  Station  (NUWES)  in  Keyport, 
Washington,  has  since  1969  been  operating  a  system  called  T-NAAS 
(Torpedo  Noise  Acquisition  and  Analysis  System).  This  system  has  later 
been  upgraded  and  today  is  called  the  Noise  Recording  System  ( NRS ) .  It 
is  used  to  measure  torpedo  radiated  noise  in  the  frequency  range  of  50 
Hz  to  20  kHz  and  the  results  have  been  used  to  establish  torpedo 
radiated  noise  baseline  levels,  to  evaluate  noise  reduction  attempts, 
and  to  derive  a  torpedo  radiated  noise-proofing  criterion  and  to  study 
parameters  affecting  torpedo  noise.  These  measurements  are  now  carried 
out  on  a  routine  basis. 

The  empirical  transmission  loss  model  that  is  currently  being  used 
for  short  range  acoustic  propagation  is  based  on  limited  empirical 
data,  but  seems  to  work  satisfactorily  under  normal  operating 
conditions.  However,  it  is  possible  under  unusually  rough  surface 
conditions  or  when  a  directional  acoustic  array  (Bottom  Mounted  Array 
(BMA))  is  used  that  this  model  would  be  inappropriate  thus  leading  to 
unknown  errors.  Both  experimental  and  theoretical  studies  have  been 
done  over  the  past  years  by  other  students  and  faculty  at  NFS  to 
improve  the  acoustic  transmission  loss  model.  Initial  experiments 
performed  in  Dabob  Bay  by  LT.  Elsen  and  Prof.  Wilson  showed  that  at 
short  to  moderate  ranges  the  contributions  to  the  signal  received  by  an 
omnidirectional  hydrophone  from  the  bottom  and  surface  reflected  paths 


were  a  significant  fraction  of  the  direct  signal  level.  They  also 
pointed  out  that  if  an  array  were  to  be  used  instead  of  an 
omnidirectional  hydrophone  some  additional  chcmges  to  the  transmission 
loss  model  should  be  made.  Theoretical  studies  presented  by  LT.  ward 
in  his  thesis  (1980)  include  the  development  of  a  computer  model  which 
calculates  the  interference  effect  in  the  acoustic  pressure  sensed  by 
an  omnidirectional  hydrophone  due  to  multipath  propagation.  This 
model,  however,  has  not  yet  been  verified  experimentally. 

The  latest  developments  in  FFT-based  dynamic  signal  analyzers  now 
maXe  it  possible  to  design  experiments  or  routine  measurements  where 
multipath  propagation  effects  on  the  test  range  may  be  studied  or 
easily  measured  and  accounted  for. 

B.  SCOPE  OP  THESIS 

In  1967,  Richard  C.  Heyser  presented  a  paper,  "Acoustical 
Measurements  by  Time  Delay  Spectrometry”  [Ref.  1],  describing  a  new 
acoustical  measurement  technique  which  would  "provide  a  solution  for 
the  conflicting  requirements  of  anechoic  spectral  measurements  in  the 
presence  of  a  reverberant  environment."  He  describes  the  technique, 
which  uses  commercially  available  instrviments,  where  the  transmitted 
signal  is  a  linear  swept  tone  (or  Linear  Frequency  Modulated  (LFM) 
pulse)  with  constant  amplitude.  This  orcing  function,  which  linearly 
relates  frequency  with  time,  provides  the  wanted  spatial  discrimination 
of  different  path  lengths  when  the  received  signal  is  analyzed  by  a 
frequency-tracking  analyzer  which  can  be  offset  in  order  to  "focus"  in 


on  the  desired  path. 


Time  delay  spectrometry  (hereafter  called  TDS)  has  since  been  used 
successfully  in  areas  such  as  investigation  of  regenerative  sound 
systems  (i.e.,  probing  of  the  sound  pressure  field  in  rooms),  free 
field  calibration  of  microphones  in  ordinary  rooms,  calibration  and 
performance  evaluation  of  miniature  ultraisonic  hydrophones,  and 
attenuation  and  velocity  meaisurements  in  tissue  as  a  function  of 
frequency.  However,  so  far  as  we  have  been  able  to  determine,  nothing 
has  been  done  in  applying  this  technique  in  the  area  of  underwater 
acoustics . 

It  is  the  scope  of  this  thesis  to  use  the  concept  of  TDS  as  the 
foundation  to  develop  a  multipath  measurement  system  which  can  be 
applied  at  the  Naval  Undersea  warfare  Engineering  Station  (NUWES)  test 
ranges  and  which  could  lead  to  a  more  accurate  computer-controlled 
model  for  transmission  loss  as  a  function  of  frequency.  In  order  to 
accomplish  this,  the  following  objectives  were  set: 

( 1 )  Identify  a  technique  for  resolving  the  direct  and  reflected 
components  of  a  sound  signal  received  by  a  hydrophone  from  a 
distant  source. 

( 2 )  Develop  the  theory  and  experimental  parameters  for  applying  this 
technique  to  measurements  on  NUWES  test  ranges. 

(3)  Develop  techniques  for  processing  the  received  signal.  This 
includes  identifying  suitable  heurdware  and  generating  necessary 
software  to  enable  processing  of  the  acquired  data. 

(4)  Verify  the  technique  and  processing  capabilities  that  have  been 
developed  by  using  selected  data  that  have  been  acquired  by 
NUWES  during  ramge  testing. 


NUWES  operates  several  ranges  for  underwater  meeisureinents .  Of 
special  interest  are  the  two  3-D  tracking  ranges,  one  at  Nanoose 
located  in  the  straits  of  Georgia,  British  Columbia,  Canada,  and  the 
other,  the  Dabob  Bay  ramge  located  off  the  Hood  Canal  in  the  state  of 
Waushington.  Even  though  the  test  data  that  have  been  used  to  develop 
the  TDS  model  was  obtained  in  Dabob  Bay,  it  will  be  shown  that 
geometrical  and  environmental  parameters  are  an  integrated  part  of  the 
model  and  hence  it  is  not  tailored  to  the  specific  conditions  of  one 
particular  range. 

Another  governing  factor  has  been  that  hardware  required  for  this 
analysis  be  compatible  with  and  the  software  be  readily  transferable  to 
existing  NUWES  amalysis  systems. 

The  above  list  provides  a  logical  description  of  the  problem  and 
will  be  the  sequence  in  which  the  subjects  are  treated  in  this  thesis 


The  "classical”  TDS  technique  as  described  by  Heyser  in  Ref.l  is 
based  on  the  use  of  a  sweeping  heterodyne  analyzer.  The  same  is  true 
for  the  more  specialized  implementation  in  the  Bruel  and  Kjaer  Time 
Delay  Spectrometry  System  Type  9950.  In  Ref. 11,  P.A.  Le'/ir,  describes  a 
calibration  and  performance  evaluation  system  in  which  a  sweeping 
heterodyne  analyzer,  Hewlett  Packard  Model  35B5A,  with  a  built  in 
frequency  offset  unit  was  used.  Common  characteristics  of  these 
systems  are  that  they  provide  excellent  measurement  results. 

However,  the  systems  are  either  too  specialized  {B«ac  Type  9950)  or  the 
analyzer  becomes  too  expensive  to  be  cost  effective  for  range  tests. 
Furthermore,  implementation  of  a  real  time  measurement  system  would 
involve  synchonization  problems  of  the  source  and  the  tracking  filter 
in  the  spectrum  analyzer. 

For  this  type  of  ajialysis,  it  is  desirable  to  record  the  data  on 
magnetic  tape  during  range  tests  and  subsequently  perform  the  analysis 
on  the  recorded  data  in  a  shorebased  facility.  This  analysis  method 
requires  a  less  specilized  (and  cheaper)  analyzer  and  the 
sychronization  problem  is  avoided  as  will  be  shown.  In  addition,  the 
test  may  be  reproduced  at  any  later  time  when  required.  This  modified 
TDS  technique  requires  that  a  Dynamic  Signal  Analyzer  be  used  to 
process  the  input  (recorded)  signal.  Dynamic  Signal  Analyzers  are 
based  on  a  high  speed  calculation  routine  which  acts  like  a  parallel 


filter  analyzer  with  hundreds  of  filters  that  are  connected  to  the 


input  signal  all  the  time.  These  analyzers  are  fast,  but  have  a 
limited  resolution  as  will  be  shown  later. 

A.  DEVELOPMENT  OF  THE  THEORY 

The  analytical  verification  of  this  modified  TDS  technique  is  based 
on  the  work  done  by  Heyser  [Ref.l]  amd  also  parts  of  the  theoretical 
verification  of  B  &  K  TDS  System  Type  9950  [Ref. 2  amd  3].  In  this 
attempt  to  bring  the  TDS  technique  to  sea,  it  will  be  necessary  to 
state  the  assumptions  that  are  inherent  in  this  linear  systems 
approach.  Also  the  nature  of  the  repetitive  LFM  (Linear  Frequency 
Modulated)  pulse  must  be  analyzed  and  finally  the  validity  of 
identifying  the  dynamic  signal  analyzer  display  with  the  acoustic 
spectrum  will  be  discussed. 

1.  Assumptions 

In  the  modified  TDS  technique,  we  consider  the  propagation  of 
small  amplitude  acoustic  signals  in  the  ocean  which  can  be  described  bu 
the  linear,  inhomogeneous  wave  equation.  Since  the  wave  equation  for 
small  eunplitude  acoustic  signals  is  linear,  we  can  represent  the  ocean 
medium  as  a  linear  filter.  Furthermore,  we  assume  that  there  is  no 
relative  motion  between  source  and  receiver  during  a  measurement  period 
and  that  the  water  column  in  homogeneous  (i.e.,  isospeed  conditions 
exist ) . 

2 .  The  Excitation  Signal 

The  excitation  signal  used  in  the  TDS  technique  has  the 
characteristic  shown  in  Fig. 2.1.  It  is  a  linear  sweep  with  constant 


The  instantaneous  frequency  is  defined  as  the  time  derivative  of  the 
phase 


i  2TI  dt  iTl  c 


(2.2) 


Note  that  the  constant  rrp/T  (or  D/2T  when  emgulair  frequency  is  used)  is 
often  referred  to  as  the  phase  deviation  constant  or  modulation  index. 
Xn  TDS,  the  value  of  the  modulation  index  is  normally  very  large.  The 
sweep  rate  S  is  the  time  derivative  of  f^^,  which  upon  differentiation 
yields : 


3  H  ^ 


(2.3) 


Prom  Eq.  (2.3)  it  is  seen  that  the  sweep  rate  S  is  constant  amd  that 
the  instantaneous  frequency  varies  linearly  with  time.  It  must  be 
noted  that  the  instaneous  frequency  f^  is  a  time  domain  description  of  the 
signal  amd  is  not  a  frequency  in  the  Fourier  sense . 

In  order  to  study  the  spectrum  of  the  excitation  signal,  we  first 
rewrite  the  signal: 


In  terms  of  the  angular  frequency,  the  phase  (p(t)  is 


In  Eq,  2.4,  each  con^Iex  ei^nential  is  separable  into  the  product  of  a 
steady^state  term  and  a  term  periodic  in  T.  Hence  it  is  sufficient  to 
expand  only  the  periodic  term  in  a  Fourier  series.  Considering  first 
only  positive  frequencies: 
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Equation  2.6  can  be  written  as 
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Conpletlng  the  square  In  the  exponent.  A,  gives 


,  1/2 

y  t|  -  nw  t  4  - 


This  can  be  put  in  the  form 
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Substituting  the  expression  for  the  exponent.  A,  gives 
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amd  chainging  the  variables  where 
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we  can  vnrite  the  Fourier  coefficient  after  multiplying  amd  dividing  by 
(2/Tr)^^^  ass 


1/2  2 
|i.j  |22j  g-KT/aO) 


T/a 


iirj 

-T/a 


a 

lx  . 
e  dx 


(2.8) 


Recalling  that  the  instamtameous  angular  frequency  is  given  by: 
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By  solving  for  t  and  svibstituting  this  into  the  expression  for  x,  Eq. 


(2.8)  can  he  put  in  the  final  form : 


where 


X 


1/2 

C"i  -  (^c 


+  ntJ^)l. 


Each  of  the  two  integrals  in  brackets  consists  of  the  Fresnel  cosine 
integral  C(u)  and  sine  integral  S(u) 
u 


a 

e‘*  dx  -  C(u)  +  i  S(u). 


(2.10) 


0 


The  f^esnel  cosine  and  sine  integrals  are  tabulated  in  Ref. 12  and  their 
behavior  is  shown  in  Fig. 2. 2.  Using  this  result  and  noting  that  the 
argument  chamges  sign  in  the  second  integral,  we  can  write  the  Fourier 
coefficients  as; 


C 

n 


2  \  ' 
Ij  ^-KT/aO)  +  C(W2>  - 
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Pig. 2. 3  shows  how  to  calculate  the  n'th  coefficient.  The  Cn's  as  given 
by  Eq.  (2.11)  are  the  side  bemd  terms  and  consist  of  the  product  of  the 
quemized  spectrum  of  a  continuous  swept  LPM  pulse  and  a  complex  Fresnel 
modifying  term  which  involves  the  finite  frequency  limits.  The  spectrum 
is  synmetrical  about  the  effective  carrier  and  even  though  there  is  an 
overshoot  at  the  band  edge,  the  ripple  rapidly  approaches  zero  as  the 
angular  frequency  is  increased.  These  ripples  are  not  distinguishable 
to  the  spectrum  analyzer  as  they  are  smoothed  by  the  appropriate  filter 
bandwidths  and  temporal  weighting  window.  This  means  that  the 
distribution  is  effectively  that  of  a  constamt  amplitude  continuum  which 
is  desired  in  TDS. 

’Aie  previous  derivation  has  revealed  several  important 
properties  of  the  forcing  function  used  in  TDS.  These  properties  may  be 
summarized  as  follows: 

-  The  spectrum  is  quite  well  behaved  for  all  modulation  indexes 
(0/2T). 

-  At  the  band  edges  or  u)^)  the  spectrum  level  is  down  6  dB 
and  falls  off  beyond. 

Expansion  of  the  negative  frequencies  yields  a  similar  symmetric 
spectrum.  This  ensures  that  the  modulation  we  are  using  is  valid 
for  TDS  even  if  the  sweep  goes  through  zero  beat. 

2 .  Response  of  the  Dynamic  Signal  Analyzer 

The  signal  received  by  a  hydrophone  when  a  lineeur  sweep  is  being 
used  to  excite  the  ocean  medium  consists  of  multiple  time-delayed 


» 


» 


» 


» 


» 


» 
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nG.2.2  THE  FRESNEL  INTEGRALS,  C(u)  AND  S(u) 


responses.  In  the  modified  TDS  technique  the  entire  signal  from  the 
hydrophone  is  sent  to  and  recorded  on  magnetic  tape  as  an  analog  signal. 
This  signal  is  then  in  turn  analyzed  by  the  Dynamic  Signal  Analyzer . 

Heyser  [Ref.l]  shows  that  (in  air)  amy  delayed  response  may  be 
selected  by  properly  adjusting  the  frequency  offset  of  the  analyzer  used 
to  process  the  transmitted  LFM  pulse.  In  the  modified  TDS  technique, 
however,  the  effect  of  a  frequency  offset  is  obtained  by  using  a  trigger 
which  is  sent  to  the  analyzer,  but  which  is  delayed  with  respect  to  the 
start  of  the  linear  swept  tone.  (Recall  that  the  sweep  rate  s  is  the 
conversion  factor. )  In  deriving  the  expression  for  the  output  of  the 
demodulator  in  a  swept  analyzer,  Heyser  [Ref.l]  shows  that  the  proper 
spectrum  is  being  displayed  and  that  the  result  is  also  valid  for 
repetitive  sweeps.  Furthermore,  the  spectrum  is  not  only  preserved  in 
amplitude  but  also  in  phase. 

B.  EXPERIMENTAL  PARAMETERS 

When  applying  the  modified  TDS  technique  to  measurements  in  the 
ocean  it  must  be  recalled  that  we  are  considering  multipath  propagation 
over  relatively  short  distances,  i.e.,  less  that  1000  yds  (approximately 
900  m).  As  a  first  approximation  the  sound  speed  profile  is  considered 
to  be  isovelocity  throughout  the  water  column  and  the  rays  will 
therefore  travel  in  straight  paths.  The  surface  and  bottom  boundaries 
are  considered  to  be  perfectly  reflecting  surfaces . 

Before  making  measurements  on  the  range,  both  geometrical  and 
environmental  parameters  need  to  be  accounted  for  in  order  to  measure 
the  effects  of  multiple  paths  between  source  and  receiver.  The  a  priori 


knowledge  about  range  geometry  and  environment  will  be  used  in  TDS  to 
determine  the  test  geometry  when  planning  measurements,  amd  subsequently 
to  identify  the  analysis  results.  As  a  tool  it  is  therefore  necessary 
to  develop  a  simple  multipath  model  to  describe  the  time  delays  caused 
by  different  path  lengths.  These  results  will  be  used  later  when 
criteria  for  the  source  signal  are  derived. 

1.  Geometrical  Considerations 

The  simple  multipath  geometry  is  shown  in  Pig. 2. 4.  Following 
the  developments  of  Albers  [Ref. 13,  p.49]  it  can  easily  be  shown  that 
the  reflected  signals  travelling  the  extra  distances  from  their  "image 
sources"  will  reach  the  receiver  at  a  time  later  than  does  the  direct 
path  sound. 

Prom  the  geometry  of  Pig. 2. 4  it  is  seen  that  the  expressions  for 
the  direct  path  (XR),  the  surface  reflected  path  (XSR),  amd  the  bottom 
reflected  path  (XBR)  aure  given  by: 


XR  -  [R^  +  (ZR  -  (2.12) 

XSR  =•  [R^  +  (ZR  +  (2.13) 

2  2  1/2 

XBR  =  [R  +  (2H  -  ZR  -  ZS )  ]  '  (2.14) 
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If  the  range  is  given  in  yards  and  the  depth  is  in  feet,  then  the  last 


term  in  Eqs.  (2.12)  to  (2.14)  must  be  divided  by  three.  The  travel  time 
for  each  path  is  now  found  by  simply  dividing  by  the  speed  of  sound,  c. 
Purthermore,  the  time  delays  between  the  reflected  and  direct  path  are 
given  by: 


Surface  reflected  time  delay:  t_s 


XSR-XR 

c 


(2.15) 


Bobtcmi  reflected  time  delay: 


t  b  -  ^®B-2£E 
~  c 


(2.16) 


and  finally  the  delay  between  the  reflected  paths: 

Inter  reflected  time  delay:  t_r  »  (2.1': 

Typical  values  for  arrival  times  and  time  delays  for  a  sound 
speed  of  4850  ft/sec  (or  1478  va/s)  are  listed  in  Tables.  I  -  X  in 
Appendix  A  on  page  77.  The  values  are  calculated  for  a  fixed  range  emd 
receiver  depth  of  300  feet  (91  m)  when  the  source  is  varied  in  depth 
from  50  to  550  feet  ( 15  -  168  m). 

2.  Environmental  Considerations 

Oceanographic  emd  acoustic  characteristics  of  both  Dabob  Bay  and 
Nanoose  3-0  test  ranges  ctre  well  documented  in  Ref. 10  and  Ref. 14.  In 
the  TDS  multipath  model  the  ultimate  goal  is  to  measure  treuismission 
loss  as  a  function  of  frequency  for  short  ranee  propagation  Also  the 
contribution  from  reflected  paths  will  be  meeisured.  Data  used  in  this 
thesis  to  develop  the  computer  controlled  multipath  model  was  taken  in 


Dabob  Bay.  The  2ureaLl  extent  of  the  range  within  Dabob  Bay  basin  is 


p 


shown  in  Pig. 2. 5  and  the  contour  of  the  basin  in  the  measurement  area 
(at  the  8000  yard  station)  is  shown  in  Fig. 2. 6. 

In  order  to  predict  arrival  times  when  plemning  an  experiment, 
information  about  the  sound  speed  profiles  are  contained  in  the  previous 
references.  For  Dabob  Bay  five  typical  profile  envelopes  have  been 
adopted  after  an  oceanographic  data  study  spemning  the  years  1971 
through  1975.  The  probability  of  occurrence  for  each  profile ‘type 
during  a  year  is  shown  in  Fig. 2. 7.  Due  to  the  lack  of  the  actual  sound 
speed  profile  for  the  measurements  performed  in  September  1982,  a  sound 
speed  profile  (envelope)  type  4  as  shown  in  Fig. 2. 11  has  been  used  in 
the  TDS  multipath  model. 

a.  Effects  of  Refraction 

AS  shown  in  Fig. 2. 9  -  2.12,  there  might  exist  a  strong 
negative  gradient  above  approximately  100  ft  (30  m).  Only  during  the 
winter  months,  ais  expected,  does  one  find  a  positive  gradient  in  the 
upper  layer  as  shown  in  Pig. 2. 8.  Another  characteristic  feature  is  the 
positive  gradient  below  300  ft  (91  m)  cis  shown  in  Fig. 2. 9,  which  is  most 
likely  to  occur  during  April  and  May  but  can  also  occur  in  late  March  or 
early  JXine. 

The  assumption  about  isovelocity  conditions  was  required  for 
the  simple  multipath  geometry  and  calculation  of  expected  time  delays  in 
particular.  In  TDS  it  is  necessary  that  a  sufficient  time  delay  exist 
between  direct  and  reflected  paths  euid  eventually  between  reflected 
paths  if  the  source  and  receiver  are  not  at  the  same  depth.  A 
queuntitative  measure  of  the  required  time  delay  will  be  derived  later. 

In  particular,  it  will  be  shown  that  under  given  conditions  the 
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the  measureinent  range  is  less  than  500  -  600  yds  (457  -  549  m).  Prom 


Pig. 2. 11  it  is  seen  that  a  worst  caise  situation  will  exist  if  both 
source  aind  receiver  are  located  at  100  ft  (30  m).  The  sound  speed 
profile  and  ray  diagram  aire  shown  in  Pig. 2. 13.  The  negative  graidient 
above  100  ft  (30  m )  is : 

g  -  (4870  -  4965)/100  -  -0.95  s“^ 


and  the  horizontal  r2mge  travelled  by  the  surface  refracted  ray  is  given 
by  [Ref. 6]: 


2AX  -  2C2  (C^g)Azl^^^  -  2C2  (4965/0.95)  *  100] 


1/2 


2045  ft. 


or  X  »  681  yds  (623  m). 

As  shown  the  surface  refracted  ray(s)  should  not  interfere  with 
measurement  of  the  direct  path  as  long  as  the  horizontal  separation 


between  source  and  receiver  is  kept  below  600  yds  (549  m). 


r  r I .  TPS  MULTIPATH  MEASUREMENT  SYSTEM 

In  the  previous  chapter,  the  LFM  pulse  with  constant  amplitude  weis 
analyzed  and  it  was  shown  theoretically  that  a  well  behaved  spectrum  is 
to  be  expected  and  that  this  spectrum  will  be  correctly  measured  by  a 
dynamic  signal  analyzer.  The  test  geometry  is  important  and  full 
advantage  can  be  taken  to  minimize  the  effect  of  unknown  parameters  by 
choosing  an  optimum  test  area  location.  The  Dabob  Bay  range  has  been 
used  in  the  previous  discussion  due  to  the  availability  of  data; 
however,  there  are  no  requirements  inherent  in  the  theory  that  prevents 
measurements  being  done  at  other  3-D  test  ranges. 

The  TDS  multipath  measurement  system  is  characterized  by  the 
following  subsystems: 

-  Range  and  test  configuration 

-  Data  processing  system 

The  range  and  test  configuration  that  is  described  here  have  been 
developed  by  Mr.  John  G.  Burwell  of  the  Acoustics  Division  at  NUWES, 
Keyport,  and  W2is  used  during  previous  TDS  measurements.  The  data 
processing  system  consists  of  both  hardware  and  software  required  to 
extract  the  desired  me2isurement  results.  The  parameters  that  need  to  be 
considered  in  order  to  resolve  the  different  paths  are  directly  related 
to  the  hardware  used  to  process  the  data.  As  a  consequence,  a 
description  of  the  experimental  setup  and  the  procedure  to  determine  the 
test  parameters  will  be  given.  Finally,  a  complete  documentation  of  the 


software  is  presented. 


A.  RANGE  AND  TEST  CXJNFIGURATION 

Data  recorded  for  this  study  were  obtained  with  the  Noise  Recording 
System  (NRS)  used  by  NUWES  in  routine  radiated  noise  measurements.  It 
consists  of  three  calibrated  omnidirectional  hydrophones  suspended  from 
a  surface  platform  and  associated  recording  instrumentation.  The  sound 
source  was  a  projector  (Type  F-56 )  suspended  from  another  platform.  The 
location  in  Dabob  Bay  wajs  near  the  8000  yards  (appr.  7300  m)  station. 

The  range  and  test  configuration  is  shown  in  Pig. 3.1.  Since  the  primary 
objective  in  these  tests  wais  to  determine  the  characteristics  of  the  new 
Bottom  Mounted  Array  system  (BMA),  conditions  were  not  optimized  for  the 
TDS  method. 

An  IRIG-B  time  code  signal  is  generated  at  the  computer  site  and 
broadcast  to  all  systems  used  in  the  test.  Prom  this  time  code  signal  a 
trigger  pulse  can  be  derived  to  trigger  the  source  signal  which  is  being 
generated  by  a  swept  ocillator.  The  period  of  this  trigger  signal 
determines  the  sweep  time  T  which  is  one  of  the  parameters  to  be 
considered  using  the  TDS  technique.  As  indicated  in  Pig. 3.1  the  IRIG-B 
signal  is  also  sent  to  the  NRS  platform  where  it  is  recorded  on  magnetic 
tape.  During  post  analysis  of  the  recorded  data  the  desired  trigger 
signal  can  now  be  extracted  (see  Chapter  III)  and  hence  a  conmon  time 
base  has  been  established. 


B.  DATA  PROCESSING  SYSTEM 
1 .  The  Hardware 

The  experimental  setup  Which  weis  used  to  implement  the  modified 
TDS  technique  is  shown  in  Pig. 3. 2.  It  consists  of  the  following 
components : 

-  HP-9836,  Computer 

-  HP-356 lA,  Dynamic  Signal  Analyzer 

-  Datum,  Time  Code  Reader  Model  9300 

-  Honeywell  5600E,  Tape  Recorder 

-  HP-7074A,  Plotter 

-  HP-2671G,  Printer 

During  a  TDS  roeaLSurement  period  analog  data  are  recorded  on 
tape.  The  complete  tape  channel  configuration  and  tape  log  is  given  in 
Appendix  C.  The  data  includes  the  individual  hydrophone  signals  and  the 
IRIG-B  (time  code)  signal.  In  Chapter  ll,  it  wais  shown  that  the  sweep 
rate,  S,  is  constant  and  that  the  instantaneous  frequency,  f^, 
varies  lineeurly  with  time.  Considering  the  distance  of  the  propagation 
paths  ( direct  or  reflected )  between  the  source  and  receiver  and  the 
speed  of  sound  in  the  medium,  it  Cein  easily  be  shown  that  a  time  delay 
is  converted  into  a  frequency  shift  which  is  proportional  to  the  sweep 
rate  euid  the  time  delay. 

The  functional  relationship  between  the  frequency  ( instantaneous ) 
of  the  transmitted  signal,  f^,  and  the  received  signal, 
fj^,  can  then  be  written  ais: 

f  ( t  -  nT )  -  f  ( t  -  nT  -  D )  where  n  -  0,1,2,3,.... 

T  R 
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.3.2  EXPERIMENTAL  SETUP 


and  D  is  the  time  delay  required  for  the  signal  to  propagate  to  the 
receiver.  If  the  relative  position  betvreen  the  source  and  receiver  is 
constant  during  a  measurement  period  then  the  time  delay  0  will  be 
constant  as  well.  The  instantaneous  received  frequency  f^^  is 
therefore  only  a  function  of  the  time  from  the  trigger  pulse.  The 
HP-3561A  features  a  post-trigger  which  will  initiate  measurement  of  the 
received  signal  at  any  desired  time  (limited  to  1023  time  records)  after 
the  trigger. 

The  appropriate  hydrophone  signal  to  be  measured  is  connected  to 
the  analyzer  input,  and  the  IRIG-B  time  code  signal  is  connected  to  the 
time  code  reader  where  a  one  pulse  per  second  ( 1  pps )  trigger  signal  is 
derived.  This  signal  is  then  applied  to  the  rear  panel  external  trigger 
input  on  the  analyzer.  The  delayed  post-trigger  is  controlled  by  the 
computer  program  and  the  accuracy  of  the  delay  increment  is 
approximately  l/lOO  ms  for  a  20  kHz  frequency  span.  (The  exact  value  of 
the  delay  increment  is  determined  by  the  frequency  span  P  and  only  a 
span  of  100  yjlz  will  result  in  an  exact  delay  time.  This  is  a  feature 
of  the  dynamic  signal  analyzer.  )  The  fact  that  meeisurements  by  the 
analyzer  now  can  be  initiated  at  a  precise  time  after  the  trigger  and 
incremented  accurately  is  the  key  to  obtaining  the  spectral  response  of 
the  ocean  medium  by  using  the  modified  TDS  technique. 

a.  The  TDS  Parameters 

Constraints  on  the  excitation  signal,  X(t),  discussed  in 
Chapter  I I  aure  determined  by  the  test  geometry  and  the  hardware  used  to 
process  the  received  signal.  In  particular  the  characteristics  of  the 
dynamic  signal  analyzer  are  important.  The  following  discusses  factors 


which  affect  the  choice  of  am  appropriate  test  geometry  (i.e.,  source 
and  receiver  depth,  maximum  horizontal  range)  based  on  a  desired 
frequency  span  P  over  which  the  response  is  to  be  measured.  A  lower 
limit  for  the  sweep  rate  S  will  be  derived  amd  the  effect  of  an 
exceedingly  high  sweep  rate  will  be  discussed. 

'The  HP-3561A  Dynamic  Signal  Analyzer  is  based  on  a  1024 
point  PFT  algorithm.  The  resulting  spectrum  (i.e.,  frequency  and 
corresponding  amplitude)  displayed  on  the  amalyzer  contains  400  points 
or  bins.  Por  a  given  frequency  span  P  the  displayed  frequency 
resolution  or  bin  separation  b  is  then  given  by: 

b  -  P  /  400  (Hz).  (3. 

Hence  if  P  -  20  kHz  then  the  bin  separation  b  =  50  Hz.  Another 
important  pairameter  regarding  the  analyzer  is  the  observation  time 
t_obs.  The  observation  time  or  time  record  length  is  the  time  during 
which  the  amalyzer  acquires  data  from  the  received  signal  and  is 
simply 

t_obs  -  400  /  P  “  1  /  b  ( sec ) .  (  3 . 

Typical  values  for  the  frequency  span  and  the  corresponding 
observation  time  using  this  TOS  technique  is  tabulated  below. 


TABLE  XI.  Frequency  Span  and  Corresponding  Observation  Time 


Or  Time  Record  Length  (HP-3561A) 


Preouencv  Span 


5  JcHz 
10  kHz 
20  kHz 
50  kHz 


Observation  Time 
80  ms 
40  ms 
20  ms 
8  ms 


Fig. 3. 3  shows  a  possible  spectrum  of  both  reflected  paths  and 
the  direct  path  a  specific  trigger  delay.  The  fact  that  the  reflected 
paths  show  up  at  a  lov«r  frequency  is  due  to  the  positive  sweep  rate  S 
(i.e.,  during  the  sweep  the  frequency  Increases).  The  time  delay 
between  the  direct  and  the  surface  reflected  path  is  t_s .  The  time  t_b 
is  the  time  delay  between  direct  and  bottom  reflected  path  and  t_r  is 
the  difference  in  arrival  time  between  the  surface  and  bottom  reflected 
paths. 

Data  were  acquired  by  the  analyzer  from  the  received  signal 
during  the  observation  time  t_obs.  Since  the  instantaneous  frequency  is 
changing  with  time  the  spectrum  of  the  signal  averaged  over  this  time 
interval  will  contain  energy  in  a  frequency  band  which  is  given 


f_spread  =  S  x  t_obs  (3.3) 

Compared  to  the  "classical"  TDS  technique  aa  introduced 
by  Heyser  [Ref.l]  and  further  developed  by  Bruel  &  Kjaer  [Ref. 2  &  3], 
this  frequency  spread  may  be  seen  ais  the  pairameter  which  determines  the 


choices  of  tracking  filter  bandwidth.  Prom  Pig. 3. 3  it  is  also  seen  that 
separation  between  the  reflected  paths  amd  the  direct  path  is  determined 
by  the  delay  in  arrival  between  the  three  paths  and  the  sweep  rate. 

Hence  we  may  write  other  sufficient,  but  not  necessary,  conditions  as: 

S  X  t_s  >  S  X  t_obs  (3.4) 

1 

S  X  t_r  >  S  X  t_obs  (3.5) 

and  S  X  t_b  >  2  x  S  x  t_obs .  (3.6) 

These  inequalities  must  be  satisfied  individually  to  resolve  all  the 
paths  Note  that  if  the  surface  and  bottom  time  delays  are  equal  we 
may  then  write 

t_delay  >  t_obs,  (3.7) 

where  t_delay  is  the  time  delay  between  the  direct  amd  reflected  paths. 

A  measurement  window  is  a  weighting  function  which  is 
applied  to  the  input  data  in  order  to  reduce  errors  caused  by  leakage 
into  adjacent  frequency  bins.  For  a  sufficiently  low  sweep  rate  such 
that  there  is  no  spreading  of  energy  over  frequency  as  measured  by  the 
analyzer,  there  is  still  a  minimum  bandwidth  that  is  defined  by  the 
frequency  span  P  and  the  measurement  window  that  is  used.  This 
bandwidth  is  given  by: 


f _band  =CxF.  (3.8) 

The  constant  C  is  normally  given  as  a  percentage  of  the  frequency 
span.  For  a  Hamning  window  ( HP-356 lA)  this  constant  is  0.00375  which 
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gives  f_band  =  75  Hz  for  a  frequency  span  of  20  kHz.  A  second  condition 
to  resolve  the  paths  is  therefore  to  ensure  that 

S  X  t_s  >  f_band  (3.9) 

and  S  X  t_r  >  f_band  (3.10) 

and  S  X  t_b  >  2  x  f_band.  (3.11) 

In  order  not  to  choose  the  sweep  rate  too  low,  we  consider  the  case  when 
t_r  =  0.  Then  the  above  set  of  inequalities  reduces  to: 

S  X  t_delay  >  f_band  =  C  x  F 

or  S>(CxP)/  t_delay.  (3.12) 

This  inequality  gives  an  expression  for  the  minimum  sweep  rate.  For  the 
HP-3561A  the  constant  C,  determining  the  equivalent  noise  bandwidth,  is 
0.00955  for  a  flat  top  window  and  0.0025  for  a  uniform  window.  The 
constraint  in  Eq.(3.12)  does  not  imply  that  high  sweep  rates  should  be 
used.  This  is  because  exceedingly  high  sweep  rates  result  in  less 
energy  being  transmitted  in  a  given  frequency  band  which  results  in  a 
lower  signal  to  noise  ratio  at  the  receiver. 

b.  Procedure  to  Determine  the  TDS  Parameters 

Based  on  the  previous  work  a  simple  procedure  may  now  be 
developed  in  order  to  fix  the  parameters  S,  F,  R,  ZS,  and  ZR  for  a  given 


measurement . 


1 )  .  First  the  frequency  span  P  over  which  the  data  is  to  be 
processed  must  be  selected.  Fixing  F  causes  the  observation  time, 
t.  obs  =  P  /  400,  to  be  fixed.  For  P  =  20  kHz  the  observation  time 
t_ obs  =  20  ms  as  given  in  Table  XI  on  page  47 . 

2 )  .  The  next  step  is  then  to  choose  a  geometry  such  that  t.  s  amd  t.  r 
are  guaranteed  to  be  greater  than  the  observation  time  t_obs.  If 
both  source  and  receiver  are  located  at  the  middle  dep>th  then  t.  r  =  0 
<t.s  =  t.b).  Fig. 3. 4  shows  the  relationship  between  time  delay  and 
range  for  a  total  depth  of  600  ft  and  sound  speed  4850  ft/sec.  The 
receiver  is  fixed  at  300  ft  and  the  source  depth  is  varied  between 
100  and  300  ft.  Mote  that  in  this  case  the  plot  of  range  vs.  time 
delay  is  symmetric  about  300  ft. 

3) .  Finally  from  the  selection  in  steps  1)  and  2),  the  minimum  sweep 
rate  must  be  determined  from  Eq .  (3.12).  Fig .3.5  shows  a  plot  to 
determine  the  minimum  sweep  rate  based  on  Eq.  (3.12)  and  C  0.00375 
(Hanning  window),  for  a  frequency  span  F  =  5,  10,  20,  and  50  kHz. 


c.  Determining  the  TDS  Parameters  —  An  Example 

Consider  a  case  where  it  is  desirable  to  cover  a  frequency 
span  of  20  kHz.  Prom  Table  11  we  get  the  corresponding  time  record 
length  or  observation  time  t_obs  =  20  ms.  In  order  to  separate  the 
direct- both- the. reflected- pethe  a  geometry  must  now  be  determined 
so  6U3  to  guarantee  that  both  t.  s  and  t  r  be  greater  than  t_  obs  and  that 
t.  b  be  greater  than  2  t^  obs . 
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If  the  receiver  is  fixed  at  300  ft  depth  then  Fig. 3.4  gives 


the 


paths  For  source  and  receiver  at  300  ft  the  maucimum  range  is 
approximately  600  yds.  At  this  range,  however,  both  reflected  paths 
have  the  same  time  delay  amd  will  contribute  to  the  reflected  signal. 

In  order  to  resolve  all  paths  Table  I  through  X  of  Appendix  A  may  be 
used  to  determine  the  delays.  From  Table  III  we  see  that  a  range  of  300 
yds  and  receiver  at  300  ft  depth,  a  sorce  depth  of  either  200  or  400  ft 
will  be  required. 

The  minimum  sweep  rate  is  found  directly  from  Fig. 3. 5.  For 
a  required  time  delay  of  20  ms  the  minimum  sweep  rate  is  approximately  4 
Hz/ms  or  equivalently  4  kHz/sec. 

2 .  The  Software  and  its  Documentation 

The  code  written  to  accomplish  the  measurements  in  the  TDS 
Multipath  Measurement  System  consists  of  2134  lines  and  is  not  included 
in  this  thesis.  The  HP  BASIC  prograiraning  language  with  extension  AP2_1 
has  been  used.  As  shown  in  Fig. 3. 2  the  program  (TDS_1,  i.e.,  version 
II^I)is  designed  to  run  on  an  HP-9836  desktop  computer.  A  modular  design 
structure  was  chosen  in  order  to  reduce  the  testing  time  of  the  code  and 
to  increase  readcUbility  and  maintaincUiility.  The  individual  modules  are 
all  subprogreuns  that  are  called  upon  «^en  required  to  perform  a  specific 
teuak.  In  addition  to  cin  easily  identifiable  structure,  comments  are 
included  in  the  subprogreun  headings  and  throughout  the  program  in  order 
to  support  readability  of  the  code.  In  particular  mnemonics  used  to 
control  the  HP-3561A  Dynamic  Signal  Analyzer  are  explained  at  the  first 


time  that  they  occur  in  a  statement. 


The  hterachical  structure  of  the  TDS_1  program  is  shown  in 
Pig. 3. 6.  Three  major  subprograms  form  the  basis  of  the  total  program 
and  represent  the  different  analysis  methods  that  can  be  used  to  process 
the  recorded  data.  The  QuicK_ Analysis  emd  Normal_Analysis  subprograms 
are  both  interactive  programs  which  require  an  operator  to  identify  the 
different  paths  (direct-,  surface-,  or  bottom- ref lected )  when  displayed 
on  the  ainalyzer.  Both  these  subprograms  aire,  however,  fundamentally 
different  with  respect  to  the  mode  of  operation  of  the  spectrum 
analyzer.  The  Auto_Analysis  subprogram  is  in  principle  an  automated 
version  of  the  Normal_Analysis  subprogram  where  path  identification  is 
done  by  the  program.  A  description  of  those  subprograms  that  support 
all  the  three  major  analysis  methods  will  first  be  discussed  followed  by 
a  description  of  each  of  the  major  subprograms  amd  their  variables, 
a.  The  structure  and  General  Purpose  Subprograms 

As  shown  in  Pig .3.6  the  only  purpose  of  the  main  program  is 
to  call  the  TDS_Henu  and  finally  to  return  the  analyzer  to  local  control 
and  terminate  the  program.  The  TDS_Henu  is  a  subprogram  that 
initializes  the  desirable  analysis  method  (i.e..  Quick,  Normal,  or 
Auto).  Before  continuing,  the  following  terminology  should  be  kept  in 
mind.  The  word  mode  refers  to  the  different  modes  of  operation  (i.e.. 
Time  Capture  Mode,  Marrow  Band  Mode,  etc.  )  of  the  Dynamic  Signal 
Analyzer.  On  the  other  hand,  the  word  method  refers  to  the  three 
different  analysis  methods  implemented  in  the  TDS_1  program. 

Regardless  of  which  analysis  method  is  being  selected  from 
the  TDS_menu,  the  Set_analyzer  subprogram  is  used  by  all  three  major 
subprograms.  Initially  the  cuialyzer  is  put  in  a  known  preset  state 


t 
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Main  program 


FIG. 3.6  THE  HIERARCHICAL  STRUCTURE  OF  THE  TDS  PROGRAM 


and  then  the  Time  Caipture  mode  which  is  used  in  the  Auto_analysis  is 


set.  If  Normal^amalysis  or  Auto_analysis  is  selected  from  the  TDS_menu 
then  the  analyzer  is  changed  to  Narrow  Band  mode  in  the  Charge_setup 
program.  After  the  HP-3561A  is  programmed  for  a  particular  analysis 
method  the  instrument  state  will  be  displayed  on  the  analyzer.  The 
Plot_display  is  at  present  only  called  in  Quick-analysis,  where  the  3-D 
magnitude  map  initially  displayed  on  the  amalyzer  after  a  mecusurement 
may  be  plotted  on  an  HP-7470  plotter. 

The  following  subprograms  are  all  used  in  connection  with 
identification  of  the  received  signal(s)  and  meaisurements  of 
corresponding  frequency  and  amplitude  values.  Quick^heading  and 
Normal_heading  are  called  in  the  Quick_analy3i3  and  Normal^amalysis 
subprograms  respectively.  In  these  subprograms  the  signal 
identification  is  entered  by  the  operator  and  the  result  is  passed  to 
the  Decision  subprogram  which  measures  frequency  and  amplitude  aund 
stores  the  values  in  the  appropriate  variables.  Set_delay  is  called  in 
both  Normal  and  Auto_analysis  and  is  used  to  increment  the  delayed 
trigger  signal.  Dump_trace  records  the  amplitude  values  corresponding 
to  the  trace  that  is  currently  displayed  on  the  analyzer  and  in  euldition 
calculates  the  average  amplitude  value.  This  average  amplitude  is 
assigned  to  the  variable  Noise_ampl  in  all  three  analysis  methods. 
However,  in  Auto_omalysis  the  Ouinp_trace  subprogram  is  also  used  to  set 
am  appropriate  threshold  determined  by  the  subprogram  Threshold.  When 
both  the  direct  and  reflected  paths  are  present  this  initial  threshold 
will  be  raised  by  Adapt_thre8hold  to  account  for  increased  signal 
levels.  The  Window_ function  is  only  used  in  Auto_analysis .  Here  a 
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measurement  window  is  moved  across  the  display  and  a  peaX  is  counted  i£ 


every  an^litude  value  within  the  window  exceeds  the  current  threshold. 
Finally  identification  of  the  peaks  are  made  based  on  trigger  delay, 
nundser  of  peaks  and  their  relative  position  on  the  display.  If  too  many 
peaks  are  counted  then  Error_me3sage  is  called  to  notify  the  operator  of 
this  event. 

The  operator's  response  to  prompts  by  the  TDS_1  program  is 
always  tested  at  the  time  the  input  occurs.  In  case  of  an  incorrect 
answer  the  subprogram  tirrong_key  is  called  to  correct  the  mistake  before 
the  next  statement  is  executed.  The  meaisurement  results  for  each 
analysis  method  are  finally  presented  in  a  tabulated  form  (Appendix  B) 
and  if  desired  a  plot  of  amplitude  vs.  frequency  for  both  the  direct  and 
reflected  paths  (see  Pig, 4. 3  and  4.4)  is  made.  During  the  analysis, 
spectral  density  levels  axe  measured  by  the  analyzer.  Both  the  print 
routine  Print_result  and  the  plotting  routine  Plot_level  intially  call 
the  subprogram  Magnitude.  In  this  subprogram  the  mecisured  spectral 
density  levels  are  converted  to  a  magnitude  level  in  dB  re  IV  before  the 
data  are  printed  or  plotted. 

b.  The  Three  Major  Subprograms 

In  version  #  1  of  the  program  both  source  and  receiver  are 
aissumed  to  be  located  at  middle  depth.  This  configuration  gives  the 
minimum  number  of  possible  states  or  combinations  of  the  signals  (paths) 
displayed  on  the  spectrum  analyzer.  The  analysis  methods  have  been 
named  to  indicate  their  functional  content  in  the  TDS_1  program. 

Quick-Analysis  is  based  upon  the  analysis  of  data  from  40 
time  records.  The  appropriate  settings  are  entered  by  the  subprogram 


Set_analyzer  and  puts  the  analyzer  primarily  in  the  Time  Capture  Mcxle. 

It  is  also  import emt  to  recognize  the  following  settings:  20  XHz 
frequency  span,  500  Hz  start  frequency  and  20%  overlap  processing.  For 
a  frequency  span  of  20  kHz  the  corresponding  time  record  length  (or 
observation  time,  see  Table  XI,  page  46)  is  20  ms.  As  a  consequence  the 
total  measurement  time  is  just  40  x  20  ms  >  800  ms  and  therefore  only 
one  sweep  is  processed. 

The  intention  with  the  Quick_analy8is  method  is  primarily  to 
provide  a  quick  overview  of  the  received  signal  as  a  function  of 
frequency  and  time,  and  to  determine  the  arrival  times  of  the  direct  and 
reflected  paths.  The  arrival  times  can  eeusily  be  determined  by 
adjusting  the  start  time  of  the  measured  data  which  is  stored  in  the 
memory  of  the  analyzer.  The  overview  of  the  received  signal  is 
presented  in  a  3-diroensional  magnitude  map  which  consists  of  60  traces. 
When  20%  overlap  processing  is  used  each  consecutive  trace  is  delayed  by 
4  ms.  (Hotel  If  the  percentage  of  overlap  is  changed  in  the 
Set_analyzer  subprogram  then  the  constant  0elay_step  must  be  chamged  in 
the  Quick^auialysis  subprogram  since  this  affects  the  time  interval 
between  successive  traces . )  Information  about  the  different  arrival 
times  is  required  as  input  parameters  for  the  Auto_auialysis  subprogram 
and  must  be  recorded  manually.  If  desired,  the  analysis  may  be 
continued  in  order  to  obtain  the  standard  meeisurement  results. 

Normal  auialysis  is  am  interactive  amalysis  method  where  the 
signal  identification  is  done  by  the  operator.  After  the  initial  setup 
is  finished,  the  analyzer  is  set  to  Narrow  Bamd  Node  in  the  Change_setup 
subprogram.  A  standard  20  kHz  span  is  set  in  the  Set_analyzer 
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subprogram.  This  corresponds  to  a  frequency  resolution  on  the  analyzer 


display  of  50  Hz.  As  was  outlined  in  Chapter  II,  the  sweep  rate  (20 
kHz/sec  in  this  experiment)  is  the  conversion  factor  between  time  delay 
amd  frequency.  The  trigger  signal  is  delayed  after  each  measurement  by 
the  value  of  the  constant  Delay_step.  This  constant  should 
theoretically  be  equal  to  or  greater  than  3  ms  which  will  give  a 
frequency  increment  for  each  delay  of:  3  ms  x  20  Hz/ms  =  60  Hz.  A  2  ms 
delay  step  will  similarly  give  a  frequency  increment  of  40  Hz  which  is 
less  than  the  resolution  of  the  spectrum  analyzer  in  Narrow  Band  Mode  at 
a  frequency  span  of  20  kHz. 

In  the  Normal_analysis  method  eight  consecutive  sweeps  are 
RMS  averaged  to  increase  the  signal  to  noise  ratio  before  a  measurement 
is  made.  The  averaging  is  also  required  to  reduce  the  effects  of  minor 
fluctuations  in  the  received  signal.  it  is  believed  that  these 
fluctuations  are  due  to  a  combination  of  scattering  from  surface 
roughness  and  small  variations  in  the  horizontal  range  between  source 
and  receiver.  Measurements  are  repeated  until  the  value  of  the  variable 
Trigger_delay  is  equal  to  or  greater  than  the  variable  Limit.  The  value 
of  Limit  is  initially  1000  (ms),  but  may  be  changed  by  the  operator  in 
order  to  analyze  just  the  time  spein  in  which  there  is  signal  present. 
After  the  analysis  is  terminated  the  result  may  be  viewed  on  the  screen 
before  a  hard  copy  is  printed.  Finally  a  plot  of  the  measured  levels 
vs.  frequency  for  both  the  direct  and  reflected  paths  is  presented. 

Auto  analysis  is  in  principle  an  automated  version  of  the 
NormaX_analysis  method  where  the  same  settings  are  used  for  the  spectrum 
analyzer.  In  this  method,  however,  the  signal  identification  and 


measurements  are  done  under  program  control.  Information  about  arrival 
times  of  the  direct  and  reflected  paths  is  required  in  order  to  simplify 
the  identification  problem.  These  inputs  must  be  obtained  from  the 
Quic>L.analysis  method  as  explained  earlier.  Since  the  analyzer  is 
operated  in  the  Narrow  Bamd  Mode,  the  same  considerations  about  display 
frequency  resolution  and  choice  of  the  constant  Delay_step  as  described 
in  the  Normal_analysis  method  are  valid  here  as  well.  Meeisurements  are 
terminated  when  the  variable  Trigger_delay  is  equal  to  or  greater  than 
the  variable  Limit  which  initially  may  be  changed  by  the  operator. 
Finally  the  measurement  results  are  presented  the  same  way  ets  described 
for  Normal. analysis . 

When  frequency  and  corresponding  amplitude  of  the  various 
paths  are  measured  during  an  analysis  period,  the  values  are  aissigned  to 
the  appropriate  variable  arrays.  In  both  Normal  and  Auto_analysis  eight 
consecutive  sweeps  are  RMS  averaged  before  the  frequency  amd  the 
amplitude  are  measured  for  a  given  value  of  the  Trigger.delay .  The 
variable  Delay_step  determines  (approximately)  the  frequency  shift  of 
the  next  measurement  which  under  stable  conditions  will  result  in  a 
monotonic  increasing  frequency.  In  case  of  an  unstable  spectrum  of  the 
received  signal  it  may  happen  that  the  frequency  meeusured  for  one  (or 
both)  8ignal(s)  at  a  given  Trigger_delay  is  less  them  the  preceeding 
one.  Also  in  the  Auto.analysis,  errors  may  occur  during  identification 
of  the  signals  due  to  incorrect  aunrival  time  supplied  by  the  operator, 
or  due  to  the  threshold  being  set  too  high  or  too  low.  In  any  of  these 
causes  the  measured  frequencies  will  not  be  monotonic  increausing.  The 
tcQnilated  output  from  Print_result  contains  all  the  meaisured  frequency 


and  amplitude  values.  In  order  to  make  a  plot  of  frequency  vs. 


amplitude  for  both  the  direct  and  reflected  paths  it  is  required  that 
the  frequency  arrays  are  monotonic  increasing.  This  requirement  is 
initially  tested  in  the  Plot_level  subprogram.  If  incorrect 
measurements  are  found  in  this  test  they  will  not  be  plotted. 


rv.  VERIFICATION  OF  THE  TECHNIQUE  AND 
PROCESSING  CAPABILITIES 

The  data  used  to  develop  the  TDS_1  program  described  in  Chapter  III 
were  obtained  in  1982  during  verification  of  the  bottom  mounted  array 
(BMA)  system  at  the  3-D  range  in  Dabob  Bay,  f^hington.  The  objective 
for  that  test  vas  to  determine  the  directional  characteristics  of  the 
BMA.  The  additional  data  obtained  during  that  test  and  made  available 
to  the  author  in  order  to  study  multipath  propagation  using  the  TDS 
technique,  is  due  to  the  foresight  of  John  G.  Burwell  who  at  that  time 
Wets  interested  in  applying  the  concept  of  Time  Delay  Spectrometry  to 
underwacer  acoustic  measurements. 

The  test  configuration  shown  in  Fig. 3.1  consists  of  two  systems. 

The  Noise  Reduction  System  (NRS)  is  the  receiving  system,  which 
consists  of  several  hydrophones  located  at  different  depths.  During 
this  experiment  the  received  signal  at  three  different  hydrophones 
located  at  240,  300,  and  340  ft  was  recorded  on  separate  channels  on  the 
tape  recorder.  The  tape  channel  configuration  amd  tape  log  is  found  in 


Appendix  C  on  page  91.  The  source  system,  am  omnidirectional  transducer 


(Naval  Research  Laboratory  Model  P  56 )  was  located  at  mid-depth  about 
300  ft.  (NRS  may  also  be  used  to  establish  transmitted  levels.  )  The 
horizontal  range  between  source  and  receiver  was  ]cept  constant  at  300 
yds  during  the  e3q>eriinent .  The  excitation  signal  was  an  LPM  pulse  at 
constant  amplitude  which  was  swept  once  each  second  covering  the 
frequency  range  from  500  Hz  to  20  kHz.  The  occillator  output  was  not 
blanked  during  the  fly  back. 

A.  THE  MODIFIED  TDS  TECHNIQUE  AND  ITS  ABILITY  TO  RESOLVE  MULTIPLE  PATHS 
Given  the  frequency  span  F  of  19.5  kHz  swept  once  every  second,  it 
follows  that  the  sweep  rate  is  s  •  19.5  Hz/ms.  Choosing  a  frequency 
span  of  20  kHz  on  the  Dynamic  Signal  Analyzer  corresponds  to  an 
observation  time  t_obs  «  20  ms  as  given  in  Table  XI  on  page  46 .  For  the 
Ceise  when  both  the  source  ^md  the  receiver  aire  located  at  mid-depth  ( 300 
ft),  the  time  delay  relative  to  the  direct  path  is  the  same  for  both  the 
surface  and  bottom  reflected  paths.  Only  separation  between  the  direct 
path  and  a  conibination  of  the  reflected  paths  cam  be  expected  for  this 
configuration  of  source  eind  receiver.  Recall  (see  Ch.III)  that  in  order 
to  resolve  these  paths  it  is  required  that  the  time  delay,  t_delay,  be 
greater  than  the  observation  time,  t_ob8.  From  Fig. 3. 4  this  requirement 
gives  a  maximum  range  of  approximately  600  yds  amd  hence  the  geometry  is 
well  within  the  limits.  The  sweep  rate  of  approximately  20  Hz/ms  is 
however  not  optimal  in  this  test.  From  Fig. 3. 5  the  minimum  sweep  rate 
is  4  Hz/ms  and  clearly  10  Hz/ms  would  be  sufficient  in  this  case.  The 
effect  of  an  exceedingly  high  sweep  rate  was  discussed  in  Chapter  III. 
However,  the  lowest  possible  sweep  rate  would  increase  the  total 


measurement  time  and  hence  increase  the  range  time  required  to  complete 
the  meaugurements . 

A  magnitude  map  obtained  from  the  Quick^analysis  method  where 
the  analyzer  is  set  to  the  Time  Capture  mode  is  shown  in  Pig .4.1  ( a )  aind 
(b)  for  two  different  start  times  relative  to  the  trigger.  Each  map 
contains  60  traces  which  represent  the  time  delay  axis  is  the  3-D  plots. 
In  Fig. 4.1  (a)  the  start  time  is  delayed  68  ms  relative  to  the  trigger 
and  shows  the  direct  path  ais  the  diagonal  signal.  The  reflected  path  is 
displayed  above  and  parallel  to  the  direct  path.  In  Pig. 4.1  (b)  the 
start  time  has  been  further  delayed  to  180  ms  showing  the  direct  path  to 
the  right  and  the  reflected  path  to  the  left.  The  conbination  of  the 
two  figures  shows  the  entire  spectrum  of  each  path  within  the  frequency 
span  of  the  excitation  signal.  In  addition  to  the  two  signals,  the 
second  harmonic  of  the  excitation  signal  which  has  travelled  by  the 
direct  path  is  also  visible.  The  presence  of  the  this  harmonic  is 
assumed  to  be  caused  by  nonlineairities  in  the  transducer  which  might 
have  been  driven  above  its  linear  range.  Another  effect  that  is  seen  at 
the  end  of  the  direct  path  is  the  increeised  spectral  level  when  the 
frequency  of  the  excitatin  signal  flys  back  to  the  start  value. 

Ideally,  the  source  should  be  turned  off  during  this  transition  when  the 
sweep  is  restarted.  The  fast  negative  sweep  in  this  case  results  in  a 
signal-to-noise  ratio  for  the  reflected  spectrum  which  is  significantly 
lower  (approximately  5  dB(V))  at  this  time  delay  than  earlier  or  later 
spectra.  Pig. 4. 2  shows  a  spectrum  of  the  received  signal  at  a  delay 
equal  to  180  ms  relative  to  the  trigger.  Clearly  a  sufficient 
signal-to-noise  ratio  is  obtained  for  both  the  direct  and  the 


reflected  spectrum  levels.  Prom  Fig. 4. 2  we  see  that  the  signal-to-noise 
ratio  is  approximately  35  dB  for  the  direct  path  (11300  Hz)  aind  16  dB 
for  the  reflected  path  (3500  Hz)  at  this  trigger  delay.  As  the  trigger 
is  further  delayed  still  better  signal-to-noise  ratios  are  obtained. 

An  example  used  to  determine  the  TDS  parameters  required  to  resolve 
all  the  paths  was  given  in  Chapter  m  (see  page  50).  By  the  use  of 
Tables  r  -  X  (Appendix  A)  it  wais  found  that  complete  separation  could  be 
obtained  at  a  range  of  300  yds  when  the  receiver  was  )cept  at  300  ft  and 
the  source  positioned  at  a  depth  of  200  or  400  ft.  For  the  test 
considered  here,  the  geometry  and  hence  the  spectral  separation  is 
marginal  and  in  addition  the  available  data  were  recorded  while  the 
source  waa  being  hoisted  towards  the  surface.  Due  to  the  constant 
motion  of  the  source  and  the  lack  of  precise  information  about  the 
source  depth  it  has  not  been  possible  to  obtain  a  well  defined  spectrum 
for  this  Cause. 

B.  THE  PROCESSING  CAPABILITIES  OF  THE  TDS  PROGRAM 

In  the  previous  section  it  was  documented  that  the  direct  amd  the 
reflected  paths  can  be  properly  resolved.  It  is  the  purpose  of  the 
TDS_1  program  to  identify  the  spectral  peaks  and  to  extract  the 
available  information  concerning  level  and  frequency  as  the  trigger 
delay  (or  meaisurement  start  time)  is  increased.  The  amalysis  results 
are  presented  in  the  same  format  for  each  analysis  method  and  consists 
of  a  plot  of  the  meaisured  levels  as  a  function  of  frequency  together 
with  a  table  of  the  measured  values  including  the  delays  at  which  the 
data  were  obtained.  In  Fig. 4. 3  and  4.4  the  meausured  levels  of  both 
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Quiok  analysis 


rig.4.3  PLOT  OF  MEASURED  LEVELS  (dD  re  1  V)  VS  FREQUENCY  (kHz) 
FOR  THE  QUICK-ANALYSIS  METHOD. 
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paths  are  plotted  as  a  function  of  frequency  for  QuicK_analysis  and 
Auto_analysis  respectively.  The  corresponding  tabulated  values  are 
given  in  Appendix  B  on  page  87 . 

At  this  stage  in  the  development  of  the  modified  TDS  technique  a 
quantitative  analysis  of  the  results  would  not  be  feeusible.  First, 
because  important  environmental  data  at  the  time  of  the  test  were  not 
available  to  the  author,  and  second,  because  such  an  analysis  should  be 
based  on  test  data  for  which  an  optimal  geometry  and  optimal  TDS 
parameters  were  used.  It  is,  however,  interesting  to  notice  that  the 
same  general  characteristics  of  the  measured  levels  as  a  function  of 
frequency  %«ere  obtained  using  either  QuicK_analysis,  Pig. 4.1,  or 
Auto_analysis,  Fig. 4. 2.  It  should  not  be  e]q>ected  that  the  two  methods 
%K>uld  give  the  same  results,  since  the  QuicK_analysis  is  based  on  the 
capture  of  40  time  records  in  less  than  one  sweep,  while  in 
Auto_analysis  8  consecutive  sweeps  are  RMS  averaged  before  measurements 
are  made.  From  Fig. 4. 2  it  is  seen  that  the  difference  in  measured 
levels  between  the  direct  amd  reflected  paths  at  the  same  frequency  is 
in  the  range  of  approximately  6  -  12  dB  (V).  However,  at  the  higher 
frequencies  vrtiere  the  signal-to-noise  ratio  is  low  due  to  the  sweep 
fly-back,  we  see  that  the  difference  in  the  measured  levels  increases. 

An  important  parameter  to  be  considered  is  the  frequency  resolution 
of  the  amalyzer.  The  value  of  this  parameter  is  inversely  proportional 
to  the  desired  frequency  spam  amd  is  also  equal  to  the  inverse  of  the 
observation  time.  By  decreasing  the  frequency  span  a  better  frequency 
resolution  will  be  obtained,  but  at  the  same  time  this  advamtage  has  to 
be  paid  for  by  an  increased  observation  time.  The  plots  shown  in 


Pig. 4. 3  and  4.4  are  not  continuous  functions  of  frequency,  but  represent 
measurements  at  discrete  trigger  delays.  Furthermore,  the  frequency 
shifts  between  consecutive  mecusurements  are  determined  by  the  value  of 
the  variable  Delay_step  which  increments  the  delayed  trigger  signal  to 
the  amalyzer.  In  the  present  version  of  the  TDS  program,  this  step  size 
is  constamt  (may  be  changed  by  the  operator)  during  a  particular 
analysis  method.  If  more  closely  paced  data  points  are  required,  the 
value  of  Oelay_step  must  be  reduced,  which  increases  the  total  analysis 
time.  However,  it  is  not  useful  to  maXe  the  delay  steps  so  small  that 
the  corresponding  frequency  shifts  are  less  than  or  comparable  to  the 
frequency  resolution  of  the  analyzer.  The  meetsurement  results  presented 
previously  for  the  Quick  and  Auto_amaly3is  methods  were  obtained  with 
constant  Oelay_step  equal  to  4  ms. 

If  we  define  the  effective  analysis  time  as  the  time  required  to 
complete  any  of  the  analysis  methods  without  including  the  time  to  plot 
or  to  print  the  data,  we  see  that  the  modified  TDS  technique  is  much 
slower  than  the  "classical"  TDS  technique  which  haua  am  effective 
amalysis  time  equal  to  the  sweep  time.  This  was  e^q^ected  since  we 
initially  preferred  to  record  the  data  on  magnetic  tape  before  the  final 
amalysis.  Using  the  present  version  of  the  TDS  program  the  effective 
analysis  time  for  the  QuicK_amalysis  is  found  to  be  15  minutes.  An 
additional  15  minutes  is  required  to  plot  the  magnitude  map  amd  finally 
5  minutes  are  required  for  the  output  plot  shown  in  Pig. 4. 3.  In 
Auto-analysis  each  data  point  is  based  on  the  RMS  average  of  8 
consecutive  sweeps.  For  this  method  the  effective  amalysis  time  was  40 
minutes.  Including  an  additional  5  minutes  to  obtain  the  output  plot 


shown  in  Pig. 4. 4,  we  get  a  total  analysis  time  of  approximately  45 
minutes  for  this  method. 

The  current  version  of  the  TDS  program  can  only  handle  the  case  where 
source  and  receiver  are  located  at  the  mid-depth.  This  test  geometry 
requires  the  smallest  number  of  possible  states  or  combinations  of  the 
received  spectra  at  any  given  trigger  delay.  If  we  let  1  denote  the 
state  where  signal  is  present  and  0  the  state  where  signal  is  not 
present,  we  may  then  tabulate  the  possible  combinations  as  follows: 

STATE  DIRECT  REFLECTED 

10  0 

2  10 

3  0  1 

4  1  1 

In  the  general  case  where  the  test  geometry  is  chosen  to  obtain 
separation  of  all  the  paths  and  in  addition  include  the  above  simple 
case  for  source  and  receiver  at  mid-depth,  the  possible  combinations  may 
be  tabulated  as  follows; 


In  the  above  table  REFLECTED  indicates  the  case  where  source  and 
receiver  are  located  at  the  middle  depth.  Implementation  of  the  general 
c8Lse  is  straightforward  in  both  of  the  interactive  amalysis  methods.  In 
order  to  account  for  the  surface  and  bottom  reflected  paths,  four 
additional  variables  would  be  required  to  store  magnitude  amd  frequency 
values.  Since  the  data  are  measured  in  the  subprogram  Decision,  the 
number  of  states  here  would  have  to  be  increaised  to  ten.  Minor  changes 
would  also  be  required  in  the  subprograms  Quick  and  Normal^heauSing  to 
include  the  cidditional  states.  In  the  Auto_analysis  method,  however, 
the  logic  used  in  the  subprogram  Window_ function  to  identify  auid  measure 
the  spectral  levels  and  corresponding  frequencies  would  have  to  be 
changed  completely.  Finally,  the  printed  format  in  Print_result  must  be 
changed  and  re-indexing  of  the  new  variables  in  Plot_level  would  have  to 
be  added  for  implementation  of  the  general  case. 


IR 


In  his  book.  Underwater  Electroacoustic  Measurements  [Ref- 5/  PP-lJ/ 
Bobber  v^ites:  "Acoustics,  pcurticularly  underwater  acoustics,  is  not  a 
very  precise  science.  The  second  significant  figure  in  a  measurement 
usually  is  in  doubt  and  the  the  third,  usually  meemingless.  Accuracies 
of  ±  1  dB  (or  about  10%  pressure  amplitude)  are  quite  acceptable  in  most 
circumstances.  At  least  part  of  the  reason  for  this  is  the  instability 
of  the  medium  and  unfavorable  environment  in  which  the  electroacoustic 
instruments  must  operate.  Except  under  the  controlled  conditions  that 
can  be  achieved  only  on  a  laboratory  bench,  water  is  not  the  simple, 
stable,  homogeneous,  quiet  and  harmless  medium  that  the  layman 
imagines . " 

All  transmission  loss  models,  designed  for  shallow  water  acoustic 
propagation,  require  statistically  averaged  input  parameters  describing 
the  sea  surface,  water  column,  and  ocean  bottom.  Even  at  NUWES  3-D  test 
ranges,  the  statistics  are  not  yet  complete.  In  order  to  describe  any 
of  the  test  ranges,  complicated  and  expensive  measurements  would  have  to 
be  carried  out.  As  an  alternative  to  the  currently  used  transmission 
loss  model(s),  it  is  believed  that  routine  measurements  can  be  carried 
out  on  the  test  ramges  in  the  future  using  the  technique  and  measurement 
system  described  in  this  thesis. 

The  multipath  measurement  system  which  hats  been  developed  is  b^LBed 
on  the  concept  of  Time  Delay  Spectrometry.  In  the  classical  TDS 
technique  a  "real  time"  response  may  be  measured  in  a  small-scale 


environment.  Implementation  of  the  TDS  concept  on  a  larger  scale  3-D 


test  range,  however,  has  required  the  development  of  a  modified  TDS 
technique  which  is  more  appropriate  for  measurements  in  the  field. 
Fundamental  requirements  are  that  the  source  and  the  receiving  systems 
remain  at  a  constant  range  aind  that  the  linear  sweep  be  triggered  at 
constant  time  intervals  during  a  measurement  period.  Given  these 
requirements  and  a  test  geometry  which  is  similar  to  the  one  used  during 
torpedo  radiated  noise  meaisurements ,  spectral  separation  of  the  paths  is 
obtained  with  a  signal-to-noise  ratio  in  the  range  of  10  to  40  dB. 
Finally,  in  the  TDS  Multipath  Measurement  System,  both  the  hardware  and 
the  software  are  compatible  with  already  existing  ainalysis  equipment  at 
the  Acoustics  Division  in  Keyport. 

The  ultimate  goal  as  outlined  in  the  introduction  is  to  determine 
the  transmission  loss  along  the  different  paths  ats  a  function  of 
frequency.  The  TDS  program  may  be  the  tool  to  accomplish  this  tasX.  A 
complete  description  of  the  receiving  system  in  terms  of  calibration 
insert  loss  and  hydrophone  sensitivity  is  readily  available  at  the 
Acoustics  Division  in  Keyport .  The  received  acoustic  level  ( RL(  f ) )  at 
the  face  the  hydrophone  referenced  to  1  MPa  and  to  the  bandwidth  of  the 
measurement  is  found  using  the  following  relationship: 

RL(f)  -  ML(f)  -  CAL(f)  -  HS(f)  +  rL(f)  -  47  -  AG 
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where : 


MIi(£)  is  measured  level  (dB  re  l  V)  in  the  meaisurement  BW 

CAXj(f)  is  measured  level  of  calibration  in  a  1  Hz  BW 

HS(f)  is  the  hydrophone  sensitivity  in  dB  re  1  /iPa 

II»(f)  is  the  calibration  insert  loss  in  dB 
AG  is  the  qain  of  the  autogain  amplifier 

During  calibration  of  the  receiving  system  a  bandlimited  ( 50  kHz )  white 
noise  signal  is  applied.  The  -47  dB  term  is  therefore  simply 
10  loq(BW).  Calculation  of  Rli(f)  is  not  included  in  the  present  version 
of  the  TDS  program.  In  future  work,  one  should  consider  including  the 
characteristics  of  both  the  receiving  system  and  the  source  in  the 
computer  model.  Without  any  changes  in  the  current  TDS  program  it  is 
further  recommended  that  a  "pre-  test  decision  package"  be  added  to  the 
program.  This  package  would  be  based  on  the  simple  formulae  for  maximum 
range  and  minimum  sweep  rate  given  in  Chapter  lit.  Extension  of  the  TDS 
program  to  include  an  arbitrary  source  and  receiver  configuration  will 
require  minor  changes  in  both  the  interactive  amalysis  methods  (Quick 
and  Normal,  analysis ) .  In  Auto,  analysis,  however,  a  complete  change  in 
the  logic  is  required.  It  is  therefore  recommended  that  future  tests 
being  conducted  on  the  range  include  data  obtained  from  an  optimal 
source  and  receiver  configuration  that  will  ensure  complete  separation 
of  all  the  paths.  When  such  data  are  available  the  simple  changes 
should  be  made  first  before  the  logic  of  the  Auto  analysis  method  is 


decided  upon. 


A  final  point  to  be  considered  is  the  transmit  response  of  the  source 
and  the  frequency  span  in  which  data  is  to  be  analyzed.  It  has  been 
shown  earlier  that  a  sweep  rate  of  20  kHz/sec  is  too  high.  When  still 
utilizing  the  1  pps  trigger  signal  derived  from  the  IRIG-B  signal,  a 
reduction  of  the  sweep  rate  to  10  kHz/sec  is  readily  available  if  the 
frequency  span  is  reduced  to  10  kHz.  Data  from  a  wider  frequency  range 
can  then  be  obtained  by  separate  measurements  each  covering  a  10  khz 


frequency  span. 


TABLE  II 


SOURCE  DEPTH  VS  TIHE  DELAY  AT  FlIED  RANGE  -  DAB08  BAY 


horizontal  range  :  200  (yds) 

NATER  DEPTH  :  600  (Ul  RECEIVER  DEPTH  :  300  «t) 


SPEED  OP  SOUND  :  4B30  (ft/sicl 


DEPTH 

(Yt)  TIR(isl 

TiSRdsI 

TlBRds) 

t,s  (is) 

t_b  (IS) 

t.r  (IS) 

50 

134.021 

143.221 

214.522 

9.20041 

80.5015 

71.3011 

100 

130.403 

148.683 

206.186 

18.2794 

75.7824 

57.503 

150 

127.519 

154.639 

198.0:5 

27.1204 

70.516 

43.3956 

200 

125. 4ia 

161.036 

190.094 

35.6183 

64.6759 

29,0577 

250 

124.14 

167.823 

182,39 

43.6827 

58.2495 

14,5668 

123.711 

174.954 

174.954 

51.2429 

51.2429 

0 

350 

124.14 

182.39 

167.823 

58.2495 

43.6827 

14.5668 

125.418 

190.094 

161.036 

64.6759 

35.6183 

29.0577 

450 

127.519 

198.035 

154.639 

70.516 

27.1204 

43.3956 

130.403 

206.186 

143.683 

75.7824 

18.2794 

57.503 

550 

134.021 

214.522 

143.221 

80.5015 

9.20041 

71.3011 

TYD 

>  Ti«(  delay  (is) 

for  the  direct  path. 

TYSR 

Tile  delay  (is) 

for  the  surface  reflected  path. 

TIBR 

:  Tile  delay  (is) 

for  the  bottoi  reflected  path. 

t.s  ' 

Tile  delay  betneen  direct  and  surface  reflected  path. 

t.b  » 

Tile  delay  betxeen  direct  and  bottoi  reflected  path. 

* 

Tile  delay  betieen  surface  and  bottoi  reflected  path. 
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TABLE  III 


SOURCE  DEPTH  VS  TINE  DELAY  AT  PIJED  RANGE  -  DABCB  BAY 

(ft) 

HORIZONTAL  RANSE  :  300 
NATER  DEPTH  :  &00 

SPEED  OF  SOUND  ;  4830 

(yds) 

(ft) 

(ft/sic) 

RECEIVER  DEPTH  : 

300 

DEPTH 

(Tt)  TJR(is) 

TZSRdsI 

TlBR(is) 

t.s  (IS) 

t.b  (is) 

t.r  (IS) 

30 

192.393 

199.103 

233.246 

100 

190.094 

203.069 

248.28 

1 !  S| 

188.127 

207.47 

241.554 

188.709 

212.281 

233.088 

9  i  il 

185.833 

217.473 

228.903 

31.6214 

B 1 9 

185.347 

223.024 

223.024 

37.4568 

37.4548 

0 

350 

183.833 

228.903 

217.473 

43.0501 

31.6214 

11.4287 

186.709 

233.088 

212.281 

48.3788 

23.5722 

22.8066 

430 

188.127 

241.334 

207.47 

53.4276 

19.3437 

34.0839 

190.094 

248.28 

203.069 

58.1867 

12.9756 

45.2111 

550 

192.393 

235.244 

199.105 

62.6326 

6.5121 

56.1403 

TXD 

=  Tile  delay  (is) 

for  the  direct  path. 

TYSR 

:  Tile  delay  (is) 

for  the  surface  reflected  path. 

TY8R 

=  Tile  delay  (is) 

for  the  bottoi  reflected  path. 

t.s  * 

Tile  delay  betieen  direct  and  surface  reflected  path. 

t.b  = 

Tile  delay  betmeen  direct  and  bottoi  reflected  path. 

* 

Tile  delay  between  surface  and  bottoi  reflected  path. 
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TABLE  IV 


SOURCE  DEPTH  VS  TIKE  DELAY  AT  FIKED  RANGE 

-  DABOB  GAY 

(ft) 

HORIZONTAL  RANGE  :  <00 
HATER  DEPTH  1  bOO 

SPEED  OF  SOUND  :  4830 

(yds) 

(ft) 

(ft/SIC) 

RECEIVER  DEPTH  : 

300 

DEPTH  (ft) 

rifl(is) 

r)(SR(is) 

TI8R(i$) 

t.s  (as) 

t.b  (IS) 

t.r  (IS) 

50 

252.735 

257.732 

303.205 

4.99484 

50.4499 

45.473 

100 

240.804 

297.345 

9.9707 

44.5295 

34.5587 

150 

2<9.3<8 

244.248 

291.773 

14.8994 

42.4243 

27.5249 

200 

2<8.28 

248.041 

284.442 

19.  7409 

38.1417 

13.4008 

250 

2<7.637 

272.173 

281.389 

24.5354 

33.7511 

9.21373 

300 

2<7.<23 

274.427 

276.427 

29.2043 

29.2043 

0 

350 

247.437 

281.387 

272.173 

33.7511 

24.5354 

9.21573 

<00 

248.28 

284.442 

248.041 

38.1417 

19.7609 

18.4008 

<50 

249.348 

291.773 

244.248 

42,4243 

14.8994 

27.5249 

500 

250.834 

297.343 

240.804 

44,5293 

9.9707 

34.3587 

550 

252.735 

303.203 

257.732 

50.4499 

4.99484 

45.473 

TJD  :  Till  diliy  (is)  for  thi  dirict  path. 

TiSR  :  Till  delay  (is)  for  the  surface  reflected  path. 

TtBR  :  Tiie  delay  (is)  for  the  hottoi  reflected  path, 

t.s  =  Tiie  delay  bitiein  direct  and  surface  reflected  path. 

t,b  =  Tiie  delay  between  direct  and  bottoi  reflected  path, 

t/  :  Tiie  delay  betiieen  surface  and  bottoi  reflected  path. 
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TABLE  V 


SOURCE  DEPTH  VS 

TIHE 

DELAI  AT  FHED  RANGE  -  DABOB  BAY 

- 

HORIZONTAL  RANGE 

:  500 

(yds) 

MATER  DEPTH  :  iOO 

(ft)  RECEIVER  DEPTH  :  300 

(ft) 

SPEED  OF  SOUND  : 

4850 

(ft/sec) 

DEPTH  (ft) 

TXRds) 

TISR(is) 

TTBRds) 

t.S  (IS) 

t.b  (IS) 

t_r  (is) 

50 

313.544 

317.586 

355.483 

4.04147 

41.9389 

37.8974 

312. 01& 

320.084 

350.516 

8.07065 

38.5001 

30.4295 

310.821 

322.894 

345.784 

12.0752 

34.9628 

22.8876 

200 

309.945 

326.008 

341.298 

16.043 

31.3129 

15.2899 

309.45 

329.413 

337.068 

19.9632 

27.6175 

7.6543 

300 

309.278 

333.103 

333.103 

23.8246 

23.8246 

0 

350 

309.45 

337.068 

329.413 

27.6175 

19.9632 

7.6543 

400 

309.965 

341.298 

326.006 

31.3329 

16.043 

15.2899 

450 

310.821 

345.784 

322.894 

34.9628 

12.0752 

22.8676 

500 

312.016 

350.514 

320.086 

38.5001 

8.07065 

30.4295 

550 

313.544 

355.483 

317.584 

41.9389 

4.04147 

37.8974 

TXD  :  Tile  delay  (is)  for  the  direct  path. 

TtSR  =  Tile  delay  (is)  for  the  surface  reflected  path. 

T)(BR  :  Tile  delay  (is)  for  the  bottoe  reflected  path, 

t.s  =  Tlie  delay  bitoeen  direct  and  surface  reflected  path, 

t.b  =  Tiie  delay  betoeen  direct  and  bottci  reflected  path. 

t,r  :  Tiie  delay  betaeen  surface  and  bottci  reflected  path. 


SOURCE  DEPTH  VS  TIHE  DELRV  AT  FIIED  RANGE  -  DABOB  BAY 


HORIZONTAL  RANGE  1  600 

(yds) 

HATER  DEPTH  :  600 

(ft) 

RECEIVER  DEPTH  :  300  (ft) 

SPEED  OF  SOUND  :  ASSO 

(ft/sec) 

DEPTH  (ft)  TlR(is) 

TiSRds) 

T)(BR(is)  t.s  las)  t,b  (is) 

t.r  (IS) 

50 

374.697 

378.085 

410.434 

3.38855 

35,7371 

32.3486 

too 

373.418 

380.188 

406.139 

6.76965 

32,7208 

2*. 9511 

ISO 

372.42 

382.556 

402.062 

10.1359 

29.6415 

19.5056 

200 

371.706 

385.186 

398.21 

13.4801 

26.5042 

13.0241 

250 

371.277 

388.073 

394.591 

16.7959 

23.3138 

6.51788 

300 

371,134 

391.209 

391.209 

20.0755 

20.0755 

0 

350 

371.277 

394.591 

388.073 

23.3138 

16.7959 

6.51788 

400 

371.706 

398.21 

385.186 

26.5042 

13.4801 

13.0241 

450 

372.42 

402.062 

332.556 

29.6415 

10.1359 

19.5056 

500 

373.418 

406,139 

330.188 

32.7203 

6.76965 

25.9511 

550 

374.697 

410.434 

378.085 

35,7371 

3.38855 

32.3486 

TID  •  Tile  delay  (is)  for  the  direct  path. 

TiSR  :  Tile  delay  (as)  for  the  surface  reflected  path. 

TiBR  s  Tile  delay  (is)  for  the  bottoi  reflected  path. 

t,s  =  Tiie  delay  between  direct  and  surface  reflected  path, 

t.b  :  Tiie  delay  between  direct  and  bottoi  reflected  path, 

t.r  :  Tiie  delay  between  surface  and  bottoi  reflected  path. 


SOURCE  DEPTH  VS  TIHE  DELAY  AT  FIIED  RANGE  -  DABQ6  SAY 
HORIZONTAL  RANGE  :  700  (yds) 

NATER  DEPTH  :  400  (It)  RECEIVER  DEPTH  :  300 

SPEED  OF  SOUND  !  4830  (ft/sec) 


DEPTH 

(ft)  TXR(is) 

TJSR(is) 

TXBR(is) 

t.s  (IS) 

30 

436.047 

438.962 

467.114 

2.91473 

434.949 

440.774 

5.82541 

434.093 

442.819 

8.72629 

200 

433.48 

445.094 

456.411 

11.6132 

230 

433.113 

447.394 

453.237 

14.4813 

300 

432.99 

430.316 

430.316 

17.3264 

330 

433.113 

433.237 

447.594 

20.1439 

400 

433.48 

456.411 

445.094 

22.9307 

450 

434.093 

459.775 

442.819 

25.6826 

300 

434,949 

463.343 

440.774 

28,3954 

550 

436.047 

467.114 

438.962 

31.0665 

TXD 

Tile  delay  (isl 

for  the  direct  path. 

TJSR 

Tile  delay  (is) 

for  the  surface  reflected  path. 

TXBR 

Tile  delay  (is) 

for  the  bottoi  reflected  path. 

t.s  = 

Tiie  delay  between  direct  and  surfa 

ce  reflected  path. 

t.b  * 

Tile  delay  betNier 

direct  and  bottoi  reflected  path. 

t/  * 

Tile  delay  betieer 

surface  and  bottoi  reflected  path. 

TftBLE  VI 11 


1  . 

SOURCE  DEPTH  VS  TIHE  DELflT  AT  FIJED  RANEE  -  DABQB  BAT 


HORIZONTAL  RANEE  I  800  (yds) 

HATER  DEPTH  I  400  (It)  RECEIVER  DEPTH  :  300  (ft) 


SPEED  OF  SOUND  :  4830  (ft/s(c) 


DEPTH  (ft) 

TIR(is) 

TISR(is) 

TlBRds) 

t.s  (as) 

t.b  (IS) 

t.r  (If 

50 

497.523 

300.08 

324.944 

2.35485 

27.441 

24.8841 

100 

494,541 

501.471 

521.413 

5.11017 

25.0318 

19.9414 

130 

495.811 

503.449 

318.445 

7.45742 

22.4339 

14.9743 

200 

495.275 

505.47 

515.444 

10.1954 

20.1893 

9.9939 

250 

494.953 

507.473 

512.473 

12.7204 

17.72 

4.99939 

300 

494.843 

310.075 

510.075 

15.2297 

15.2297 

0 

350 

494.953 

512.473 

507.473 

17.72 

12.7204 

4.99939 

400 

495.275 

515.444 

505.47 

20.1893 

10.1934 

9.9939 

450 

495.811 

518.445 

303.449 

22.4339 

7.45742 

14.9743 

500 

494.541 

521.413 

501.471 

25.0518 

5.11017 

19.9414 

550 

497,523 

324.944 

500.08 

27.441 

2.55485 

24.6841 

TXD  :  Tilt  delay  (is)  for  the  direct  path. 

TKSR  :  Tiie  delay  (is)  for  the  surface  reflected  path. 

TIER  s  Tile  delay  (is)  for  the  bottoi  reflected  path, 

t.s  :  Tile  delay  betiieen  direct  and  surface  reflected  path, 

t.b  :  Tiie  delay  betNcen  direct  and  bottoi  reflected  path. 

t,r  :  Tiie  delay  betueen  surface  and  bottoi  reflected  path. 


I 


TABLE  II 


SOURCE  DEPTH  VS  TIKE> DELAY  AT  FIIED  RANGE  -  DAGOB  BAY 


HORIZONTAL  RANSE  :  900  (yds) 
MATER  DEPTH  ;  600  (ft) 


RECEIVER  DEPTH  :  TOO  (ft) 


SPEED  OF  SOUND  :  4820  (ft/sic) 


DEPTH  (ft)  TlR(is) 


TISR(is) 


TlBR(is) 


t  s  (IS) 


t  b  (IS) 


t  r  (IS) 


SOURCE  DEPTH  VS  TIME  DELAY  AT  FIIED  RANSE  •  DABOB  BAY 


HORIZONTAL  RANSE  1  1000 

(yds) 

HATER  DEPTH  :  600 

(ft) 

RECEIVER  DEPTH  :  300  (ft) 

SPEED  OF  SOUND  :  4830 

(ft/sec) 

DEPTH  (ft)  TIR(es) 

TlSR(is) 

TYBRds)  t.s  (is)  t.b  (is) 

t.r  (as) 

SO 

620.701 

622.732 

642.906 

2.03096 

22.2043 

20.1536 

100 

619.93 

624.031 

640.172 

4.10132 

20.2426 

16.1412 

ISO 

619,329 

623.477 

637.394 

6.1474 

18.2645 

12.1171 

200 

618.9 

627.039 

633.172 

8.18842 

16.2716 

8.08319 

2S0 

613.643 

628.866 

632.909 

10.223 

14.2661 

4.04309 

300 

618.537 

630.807 

630.807 

12,2498 

12.2498 

0 

3S0 

618.643 

632.909 

628.866 

14.2661 

10.223 

4.04309 

400 

618.9 

633.172 

627.089 

16.2716 

8.18842 

8.08319 

4S0 

619,329 

637.394 

623.477 

18.2645 

6.1474 

12.1171 

SOO 

619.93 

640.172 

624.031 

20.2426 

4.10132 

16.1412 

5S0 

620.701 

642.906 

622.732 

22.2045 

2.05096 

20.1536 

TXO  3  Tiic  dslay  (is)  for  thi  diroct  pith. 

TXSR  3  Tits  delay  (as)  for  the  surface  reflected  path, 

TIER  3  Tile  delay  (is)  for  the  bottoe  reflected  path, 

t.s  3  Tile  delay  betueen  direct  and  surface  reflected  path, 

t.b  3  Tile  delay  betieen  direct  and  bottoi  reflected  path, 

t.r  3  Tile  delay  betieen  surface  and  bottoi  reflected  path. 


APPENDIX  B 


TABUIATED  MEASOREMENT  RESULTS 


RCCOROeO  KXCNmiDC  LEVEIJS  (d6  C«  IV) 


DIRECT  M]U3f  DIRECT  PRQ(  HZ  >  REEL  tOiOH  REPL  PRQ(  BZ  )  NOISE  ttAOl 


18  63 

17950 

-'29.00 

10050 

-55 . 98 

17.74 

18300 

-29.11 

10500 

-55 . 89 

16.  as 

18700 

-26.86 

llOOO 

-51.00 

16.13 

18950 

-27.10 

11400 

-46.46 

20.60 

18950 

-26.92 

11750 

-47 . 15 

36.05 

19050 

-26.88 

12050 

-S3 . 85 

19.75 

0 

-26  28 

12450 

-59.30 

19  75 

0 

-26.82 

12800 

-59.43 

19  75 

0 

-27.39 

13200 

-59.90 

19. 75 

0 

-28.92 

13700 

-59.67 

19.75 

0 

-27  31 

14150 

-50 . 90 

19.75 

0 

-26.33 

14350 

-58 . 79 

19.75 

0 

-25.68 

14700 

-58.28 

19.75 

0 

-26.61 

15250 

-59 . 51 

19.75 

0 

-27  19 

15500 

-58 . 47 

19. 7S 

0 

-27  93 

15050 

-57  76 

19  75 

0 

-27  77 

16400 

-58.07 

19  75 

0 

-28.53 

16800 

-58  28 

19.  /5 

0 

-28.72 

17050 

-50.29 

19  75 

0 

-30  07 

17400 

-58.27 

19.  75 

0 

-30.81 

17900 

-57.92 

19.75 

0 

-31 . 38 

18150 

-57 . 26 

19.75 

0 

-31.64 

18500 

-56.56 

19.75 

0 

-31  40 

18900 

-54.70 

19  75 

0 

-36  55 

19050 

-55  53 

19  75 

0 

19.75 

0 

-59.43 

ANALYIS  METBOOi  Atm)  ANALYSIS 


DATE  OP  ANALYSIS!  11/17/84 


RECOROeO  MAOIITUDE  LEVELS  (dB  r«  IV) 

DELAY(KS)  DIRECT  KAGN  DIRECT  PRQ(  HZ )  REEL  HA(3I  REEL  PRQ(aZ)  NOISE  MAGN 


RKCOROED 

OeiAY(lO)  OERSCT  lOUX 

MlU3frrUD£  LEVELS  (dB  r«  IV) 

DIRECT  PR0<H2)  rRQ(  HZ  > 

NOISE  Mvai 

244 

-19.42 

17500 

-27  59 

9650 

-45.30 

248 

'18. 31 

17950 

-26.90 

10250 

-45.34 

252 

-18.07 

18250 

-27  06 

10350 

-44.49 

256 

-16.97 

16650 

-25.08 

10950 

-40.60 

260 

19.75 

0 

19.75 

0 

19.75 

264 

19.75 

0 

-24.85 

11750 

-37.12 

268 

19.75 

0 

-26. 44 

12100 

-37  50 

272 

19.75 

0 

-26.33 

12250 

-43.93 

276 

19. 75 

0 

19.75 

0 

19.75 

280 

19.75 

0 

-26.34 

13150 

-47 . 90 

284 

19.75 

0 

19.75 

0 

19.75 

288 

19.75 

0 

19  75 

0 

19.75 

292 

19.75 

0 

19  75 

0 

19.75 

296 

19.75 

0 

-27  66 

14750 

-48.48 

300 

19.75 

0 

-27.33 

15000 

-46.72 

304 

19.75 

0 

19.75 

0 

19.75 

308 

19.75 

0 

-28.01 

16000 

-48.39 

312 

19,75 

0 

-29.61 

16250 

-47.58 

316 

19.75 

0 

-29  64 

16450 

-47  12 

320 

19.75 

0 

-29.33 

16950 

-46.91 

324 

19.75 

0 

-34.56 

17500 

-47.62 

328 

19.75 

0 

-32.49 

17950 

-47.44 

332 

19.75 

0 

-32.38 

18200 

-47.43 

336 

19.75 

0 

-34.02 

18500 

-47.33 

340 

19.75 

0 

19.75 

0 

-53.30 

344 

19.75 

0 

19.75 

0 

-54.15 

90 


APPENDIX  C 


TAPE  CHANNEL  CONFIGURATION  AND  TAPE  LOG 


j 

CODE 

CONTENT 

1-7 

Not  used 

8 

(D) 

Reference 

9 

(D) 

Gain  code 

10 

(D) 

Top  hydrophone 

11 

(D) 

Middle  hydrophone 

12 

(D) 

Bottom  hydrophone 

13 

(D) 

IRIG-B 

14 

(D) 

Voice  coimnents 

TAPE  LOG 

AUTOGAIN 

A. 

Calibration 

Signal 

0  dB 

B. 

Ambient  noise 

70  dB 

C. 

S%^pt  tone. 

cons tern  t  depth 

10  dB 

D. 

Swept  tone. 

varied 

depth 

10  dB 

Receiver  depth:  Top  240  ft 

Middle  300  ft 
Bottom  340  ft 

Scource  depth:  (C)  300  ft  fixed 

(D)  Varied  from  300  ft  to  surface 
Source/Receiver  horizontal  separation  nominally  300  yds. 
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